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D , Abstract 



We present the analysis of the complete glueball spectrum for the AdS'^ black hole 
supergravity dual of QCD4 in strong coupling limit: g'^N 00. The bosonic fields 
in the supergravity multiplet lead to 6 independent wave equations contributing to 
gluebaU states with J^'^ = 2++, 1+^,1 , 0++ and 0"+. Wc study the spectral 
splitting and degeneracy patterns for both QCD4 and QCD3. Despite the expected 
limitations of a leading order strong coupling approximation, the pattern of spins, 
parities and mass inequalities bare a striking resemblance to the known QCD4 glueball 
spectrum as determined by lattice simulations at weak coupling. 
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1 Introduction 



The Maldacena duality conjecture Q and its further extensions ^ state that there is 
an exact equivalence between large N conformal field theories in d-dimensions and string 
theory in AdS^^"*" x M. Subsequently Witten |^ suggested how to break explicitly the 
conformal (and SUSY) symmetries to arrive at a dual gravity description for SU{N) 
quarkless QCD^. Thus we may have at last a definite proposal for the long sought "QCD 
string". As anticipated by 't Hooft the dual correspondence for the large expansion for 
SU{N) Yang-Mills theory is the perturbative expansion of string theory. Still this theory 
is difficult to formulate, let alone solve. At present explicit calculations also require taking 
the strong coupling limit, g'^N oo, where the string tension goes to infinity {a' 0) 
and the dual theory is classical gravity. 

In this paper we complete the study of the glueball spectrum for the strong coupling dual 
description of QCD^. For comparison the analogous spectrum calculation is presented 
for QCD3, which shows a very similar pattern, which in qualitative terms can be traced 
to the underlying flat space T duality between type IIA and IIB string theories which in 
turn are the AdS duals to QCD^ and QCD^ respectively. The goal is to learn more about 
the AdS/Yang-Mills correspondence by comparing the AdS strong coupling spectrum with 
the rather well determined glueball spectrum in lattice gauge theory. Of course, the 
strong coupling expansion at best can provide a rough guide to the underlying physics. 
Nonetheless the correspondence to the continuum (i.e. weak coupling) limit of the lattice 
spectrum is surprisingly good. This comparison may prove useful to provide support for 
the conjectured Maldacena duality and to give specific information on the major strong 
coupling artifacts that must be removed as one approaches universality at the ultraviolet 
fixed point. 

In addition to new spectral calculations, we summarize the earlier work by many authors [^, 
0, |, |, 0, 0, m. In particular, we extend our earlier paper on the tensor glueball for 
QCD3 on an AdS^ black hole background to the physically relevant case of QCD4. For 
QCD3, we found that the tensor spectrum (2++) was degenerate |jl^ with dilaton(0"'"^) 
and the axion (0"' ). However, the mass gap set by a lower scalar glueball obeys the 
inequality, 

m(0++) < m(2++) . (1) 

This scenario is repeated for QCD4. The dilaton mode (0"''^) coupling to Tr[F'^] remains 
degenerate with the tensor (2^"*"), but the axion (0"' ) is heavier, consistent with lattice 
results. Due to the trace anomaly, the lowest mass scalar (0"*"^) couples to the energy 
density Tqo and again it obeys this inequality above, Eq. (|l]). As in our earlier work. 
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the goal is to see the details of the spin structure of the lowest glueball states, which 
we believe is most sensitive to the underlying gauge theory. We find that our analysis 
combined with all the earlier results can give a systematic and complete strong coupling 
glueball spectrum. 

It is useful to end this introduction with a rough overview of our results and the organi- 
zation of the paper. 

The geometrical construction for QCD4 is roughly as follows. One starts with 11 di- 
mensional M theory on AdS''' x S'^. The seven dimensional AdS''^ we take to have a 
radial co-ordinate r and Euclidean space-time co-ordinates xi, 2:2, 2:3, 3:4, rcs and xu. The 
"eleventh" dimension is taken to be compact, reducing the theory to type IIA string the- 
ory. Matter at the center of this space (r = 0) consists of N D4-branes (or NS 5-branes 
wrapping in the 11th coordinate) with world volume co-ordinates xi ■ ■ ■ x^. The 5-d 
Yang- Mills CFT "living on" the brane is dimensionally reduced to QCD^i by raising the 
"temperature", /?~^, in a direction X5 = r, parallel to the brane. The new metric is an 
AdS'^ black hole with xn compact . Compact directions on S*^ will be denoted by x^, 
a = 7,8,9,10. 

The strong coupling glueball calculation consists of finding the normal modes for the 
bosonic components of the supergraviton multiplet in the AdS^ x black hole back- 
ground. We are only interested in excitations that lie in the superselection sector for 
QCD4. So we can ignore modes for all non-trivial harmonics in that carry a non-zero 
R charge and all Kaluza-Klein (KK) modes in the two circles with U(l) KK charges. 
Imposing these restriction and exploiting symmetries of the background metric reduce the 
problem to six independent wave equations, referred to as T4, V4, A'4, M4 and L4 in 
the text. In Fig. || (left side), we plot the low mass states for each equation, labeling the 
quantum numbers for each level. We can identify the modes with the bosonic components 
of the zero mass sector of type IIA string theory: the graviton (G), the dilaton (</>), the 
NS-NS 2-form (B) and RR 1- & 3- forms (C'(i), C'(3)). The spin degeneracy of the spectrum 
is due to a spurious 0(4) symmetry of the strong coupling approximation that combines 
the 11th and three spatial co-ordinates. However, as we explain in the text, all extra 
states not observed in the lattice data for the glueball spectrum, actually carry a discrete 
"r-parity" that places them, like KK modes, outside the QCD superselection sector. 

We have included in Fig. |^ (right side) the full spectrum for QCD^. The entire analysis is 
very similar. Starting with type IIB strings in AdS^ x which is dual to = 4 SUSY 
Yang-Mills, one introduces a compact thermal circle forming an AdS^ black hole. Again 
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Figure 1: The AdS glueball spectrum for QCD/^ (left) and QCD3 (right) with mass 
eigenvalues, m„, plotted for each of the six equations labeled on the horizontal axis. 

there are six independent wave equations, labeled by 53, T3, V3, A's, M3 and L3. They 
correspond to fluctuations for type IIB fields: the graviton (G), the dilaton {(p), the NS- 
NS 2- form (B) as before and RR 0- &: 2-forms (C(o), C(2))- For both cases we also include 
volume fluctuations in the compact sphere and for type IIA and IIB respectively. 

In Sec. 2, we give general arguments for the spin and degeneracy of glueball states for 
QCD^i followed in Sec. 3 by the analysis for QCD^. For each, we give the resultant six wave 
equations and numerical values for the first ten levels. (Derivations for these equations are 
explained further in Appendix A.) For all but the lowest eigenvalue, the glueball masses, 
rrin, are well approximated by the WKB expansions [|^, |l^: = -\- 5 n + 7). 

For the first level (n = 0), we provide a simple but reasonably tight variational upper 
bound [Q. (See Appendix B for details.) 

In Sec. 4, we give the parity and charge conjugation quantum numbers of the glueball 



4 



states using the Born-Infeld action to determine the quantum numbers of the couphngs 
between gravity fields and gauge fields. The striking similarity between the QCD4 and the 
QCD3 spectra (see Fig. |l|) can be understood qualitatively in terms of T-duality, which 
relates -D4 branes in IIA to D3 branes in IIB. 

In Sec. 5, we compare the AdS strong coupling spectrum with the well determined levels 
from lattice QCD4 and remark on the relationship to the constituent gluon picture. No 
extra states are present in the AdS spectrum that couple to QCD operators, although 
the absence of the low mass 2 state is noted. We also show how the strong coupling 
expansion for the Pomeron intercept may be used to provide an estimate of the coupling 
at the crossover between the strong and weak coupling regimes. 



2 Glueball Spectrum for QCD4 



To approach QCD^ one begins with M theory on AdS'^ x S^. We compactify the "eleventh' 
dimension (on a circle of radius Ri) to reduce the theory to type IIA string theory and 
then following the suggestion of Witten raise the "temperature", with a second 

compact radius R2 in a direction r, with (3 = 27ri?2- On the second "thermal" circle, the 
fermionic modes have anti-periodic boundary conditions breaking conformal and all SUSY 
symmetries. This lifts the fermionic masses and also the scalar masses, through quantum 
corrections. The 't Hooft coupling is g^N = 27rgsNls/R2, in terms of the closed string 
coupling, gs and the string length. Is ■ Therefore, in the scaling limit, g'^N 0, if all goes 
as conjectured, there should be a fixed point mapping type IIA string theory onto SU{N) 
pure Yang-Mills theory. 

We consider the strong coupling limit at large N, where the string theory becomes classical 
gravity in the AdS"^ black hole metric, 

ds' = ir' - \)dr' + ^ dx^ + (r' - + \dnl , (2) 

^ 1=1,2,3,4,11 ^ 

with radius of curvature, R^^g = SirggNlg. We have removed all dimensionful parameters 
in the metric by a normalization setting RAdS = 1 /? = 27r/3. 
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2.1 Spin and Degeneracy of Glueball States 



In M theory the supergraviton is a single multiplet in 11-d with two bosonic fields - a 
graviton, Gmn, and a 3-form field, Amnl, as designated in Table 0. After restricting 
all indices and co-ordinate dependence to AdS"^ , we have a graviton, G^y, a dilaton 
and an NS-NS tensor field B^y. In addition there are two RR fields, a one- form G^ 
and a three-form C^y\. Furthermore, we will also consider the scalar modes coming 
from "volume" fluctuations for S*^. The relationship between M theory and IIA string 
theory nomenclature, after restricting to the AdS"^ subspace, is presented in Table |l[ The 
table gives the J^*-" quantum numbers for all glueball states. The pattern of degeneracy 
(explained below) is indicated by the rows ending with the lowest eigenvalue for each of 
the six wave equations, Eq. @: T4,V4, etc. 

The task is to find all the quadratic fluctuations in the AdS''^ black hole background that 
might survive for QCD4 in the scaling (weak coupling) limit, ignoring any Kaluza-Klein 
mode in compact manifolds (compactified for xu, for r and the spheres 5^). They 
are charge states in their own superselection sector that are clearly absent in the putative 
target theory. Additional "spurious" states will be discussed in Sec. 4 where we treat 
discrete symmetries. 
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Table 1: IIA Classification for QCD4. Subscripts to J^^ designate Pt- = —I . 



To count the number of independent fluctuations for a field of given spin, we adopt the 
following method. We imagine harmonic plane waves propagating in the AdS radial di- 
rection, r, with Euclidean time, X4. For example metric fluctuation, 

G^i/ = g^i/ + /i^!y(x), (3) 
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in the background, g^,y, are taken to have the form, /i^i/(r, X4). There is no dependence 
on the spatial co-ordinates, Xi = {xi, X2, X3, xu) and the compactified "temperature" 
direction r. 



2.1.1 Metric fluctuations 



The four dimensional field theory lives on the hypersurface co-ordinates, xi, X2, X3, 2:4, 
(with X5 = r and xn as the two compactified coordinates.) A graviton has two polarization 
indices. If we were in flat space time, we could go to a gauge where these indices took 
values only among (xi, X2, X3, xn, r) and not from the set (r, X4). The polarization tensor 
should also be traceless. This leaves (5 x 6)/2 — 1 = 14 independent components. In 
the AdS space time, we can count the number of graviton modes the same way, though 
the actual modes that we construct will have this form of polarization only at r — > 00; 
for finite r, other components of the polarization will be constrained to acquire nonzero 
values (p^ j. 



We identify the spin content of these 14 components in two steps. First note that the 
background metric is 50(4) symmetric. (It is flat in the first four of these directions, 
511 = 522 = 533 = 511,11 = r"^, while it is "warped" in the r direction, g-rr = r'^ — l/r'^). The 
system therefore has SO (4) symmetry leading to three distinct equations corresponding 
to 9, 4 and 1 dimensional irreducible representations under 5*0 (4). In Table |l], these are 
denoted by V4, and 84 respectively. 



These representations lead us to a degenerate spectrum of spins under the physical SO (3) 
rotations in xi,X2,X3, which we list below: 



• 9-dimensional representation breaks into 5 -|- 3 -|- 1 under 5*0(3), 

Gij : hij - \5ijhkk / spin-2, 

d : hi^u = hu^i / spin-1, 

(/) : hn^n = -3hn = -3/i22 = -3/i33 / ^ spin-0, (4) 

• 4-dimensional representation breaks into 3-1-1 under 50(3), 

: h^-i = hij- 7^ — > spin-1, 

Cr ■ /it,ii = hor / spiu-O, (5) 

• singlet under 50(3), 

Gtt : hrr = -4/iii = -4/i22 = -4/133 = -4/111,11 7^ ^ spin-0, (6) 
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where i,j, k = 1,2, 3. 

In addition there is a scalar field, G^, coming from the metric on the sphere (with 
m\^g = 72), which is referred to as L4 in Table ||. 



2.1.2 Three-form fields 



The behavior of the three- form field is discussed briefly in the Appendix, but we recall some 
of the main features here. The wave equation has a topological mass term which results 
in the equation being factorized into two first-order equations, yielding upon iteration two 
second-order equations. One solution is a massless 3- form field, which has solutions that 
are pure gauge for the case when there is no dependence on the sphere 5^ , and is thus to be 
ignored. The other field gives a second-order equation with m^^g = 36, but the fact that 
we have a first-order equation as the primary equation reduces the degrees of freedom 
effectively to those of a massless 3-form field. Let the propagation directions be again 
[x4,,r). If we consider the component with indices j4i23 then we get a specific nonzero 
value also for the components ^rr.ii and ^4t,ii. We can count the independent degrees of 
freedom by looking only at components that do not have the propagation directions X4,r 
among the indices. Thus we get the fields listed in Table |l[ Reducing the 50(4) states 
under rotations in xi,X2,x-i yields: 



• 4-dimensional representation breaks into 3 + 1 under 50(3), 

Bij : Aij^ii 7^ ^ spin-1, 

C123 : Aijk 7^ ^ spin-0, (7) 

• 6-dimensional representation into 3 + 3 under 5*0(3), 

Sir : Air^u 7^ ^ spin-1 , 

Cijr ■■ Aijr 7^ ^ spin-1 . (8) 

The field equations for these states have amplitudes A'4 and M4 as listed in Table |l[ 



2.2 Wave Equations and QCD4 Glueball Spectra 



The wave equations for the metric fluctuations for QCD4 have been obtained in Ref. by 
analyzing the linearized Einstein equations about the AdS^ x black hole background 
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which leads to three independent equations, T4, V4 and 54 [10, 12]. Here we complete the 
spectral analysis giving the numerical values for all glueball masses. Fluctuations A'4, M4 
and L4 can be found similarly, leading again to all together six independent equations for 
QCD/i^ expressed in a manifestly hermitian form: 
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(9) 



We shall provide a more detailed discussion on how these equations can be obtained in 
Appendix A, while concentrating here on establishing our normalization convention. 



Consider first metric perturbations of the form 

V = ^A^y^'"' > (10) 

with all other fields set to zero. We shall further fix gauge to h^^ = 0, and from the 
linearized Einstein's equation, we determine the discrete spectrum with k^^ = im. Because 
of the 50(4) symmetry in xi,X2,X3,xii, the system is highly degenerate. Three distinct 
equations for various perturbations can be obtained by the following procedure. 



Tensor: There are five independent perturbations which form the spin-2 representations 
of 50(3): 

h,, = qi,r^T^{r)e~^-\ (11) 
where i,j = 1,2,3 and qij is an arbitrary constant traceless-symmetric 3x3 matrix. 



Vector: Consider perturbations: 

/i^r = /iri = 9iv/r6 - lV4(r)e-'""^ (12) 

where i = 1,2,3 and qi is an arbitrary constant 3-vector. Both equations for T4 and V4 
have also been obtained in Ref. Q by considering the corresponding degenerate scalar 
modes. 
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Scalar: The analogous scalar perturbation is 

hrr = (r^ - r-^)Si{r) e"'"^* 



(13) 



Three-form and Volume Scalar: Next we turn to 3-form fields. It is sufficient to 
consider 

AijAi = Bij = e»jfc gfcr^A^4 (r)e~"""'" , 

Air,n=Bir = qiV^^M,{r)e-"'^\ (14) 
where qi is again an arbitrary constant 3- vector. 

Note that the metric in the direction xu is the same as that in the directions i,j,k, so 
the above functions M4 will also give the solutions for Cijk and Cijr respectively. 
(Fluctuations for Bij have been considered previously in Ref. |^]. ) 

Lastly, for the volume perturbation, we consider Q 

/i^ = L4(r)e-"^"^ 



Equation: 
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Table 2: The mass spectrum, m^, for QCD^ Glueballs 



To calculate the discrete spectrum for each of these equations, one must apply the correct 
boundary conditions at r = 1 and r = co. This issue has been discussed in several earlier 
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papers |T^. The boundary conditions are found by solving the indicial equation. In ah 
cases the appropriate boundary condition |^] at r = 1 is the one without the logarithmic 
singularity. At r = oo the least singular boundary is required to have a normalizable 
eigenstate. (See Appendix B for a listing of all boundary conditions.) Matching boundary 
conditions from r = 1 and r = oo results in a discrete set of eigenvalues m^, where n is 
the number of zeros in the wave function inside the interval r S (l,oo). We solved the 
eigenvalue equations by the shooting method, integrating from ri ~ 1 to large r^o — oo. 
The resultant spectrum is given in Table |. 

To further check our results, we have compared our numerical masses to the WKB ap- 
proximations, 

ml ~ ixl{n'^ + (5 n + 7) + 0(i) , (15) 

n 

where ^\ = 367r(r(2/3)/r(l/6))^. For each equation, individual integer constant, 5, was 
determined analytical and constant 7 was fit to the numerical data in Table ^ (See 
Appendix B). The fits to the WKB formula are accurate to better than 0.1 % for all but 
the lowest (n = 0) mode. For each lowest mode, we have also carried out an independent 
simple variational estimate. We note that in case, our numerically calculated value for uiq 
is always close and respects the variational estimate as an upper bound. (See Table ^) 



3 Glueball Spectrum for QCD^ 



For QCD^ the construction of the supergravity dual begins with Maldacena's conjecture 
for type IIB string theory in a AdS^ x background metric. Here the dual theory is 
conjectured to be the conformal field theory for = 4 SUSY SU(N) Yang-Mills in 4-d. 
The AdS curvature is induced by N units of charge on N coincident D3 branes giving 
rise to a constant volume 5-form F(5) in the product manifold. The co-ordinates in the 
AdS^ space, we label by one "radial" co-ordinate r and four space-time co-ordinates, x^, 
/i = 1, 2, 3, 4, parallel to the D3 branes. The remaining five co-ordinates in are labeled 
by Xa, a = 6,7, 8, 9, 10. 

Following the suggestion of Witten for obtaining a supergravity dual to QCD3, we break 
conformal and SUSY symmetries, by introducing a compact "thermal" co-ordinate X4 = r 
with anti-periodic boundary on for the fermionic modes. The resultant metric at high 
temperature is an AdS^ black hole, 

ds^ = (r'-}^)dT^+r' ^ dxl + {r'-\)-^dr' + dnl, (16) 

^ j=l,2,3 ^ 
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with radius of curvature, R\^g = ^irggNl^ and 3-d Yang- Mills coupling, g^N = 2gsN/R 
in terms of the string coupling Qs, string length Ig and compact circumference f3 = 
2ttR. We have removed all dimensionful parameters from the metric by adopting a simple 
normalization with RacLS = 1 and /3 = vr, for the circumference of the thermal circle. At 
high temperature (or equivalently low energies), IIB string theory in this background is 
conjectured to equivalent to QCD3. 

3.1 Spin and Degeneracy of Glueball States 

Type IIB string theory at low energy has a supergravity multiplet with several zero mass 
bosonic fields: a graviton, Gf^u, a dilaton cf), an axion (or zero form RR field) C and two 
tensors, the NS-NS and RR fields Bfj^i, and Gfj,^ respectively. In addition there is the 4-form 
RR field C(4) that is constrained to have a self-dual field strength, F(5) = (iC(4) . 

Now the task is to find all the quadratic fluctuations in the above background metric whose 
eigen-modes correspond to the discrete glueball spectra for QCD3 at strong coupling. We 
are only interested in excitations that lie in the superselection sector for QCD^. Thus for 
example we can ignore all non-trivial harmonic in that carry a non-zero R charge and 
all Kaluza-Klein (KK) modes in the thermal circle with a U(l) KK charge. The result 
of these considerations, discussed in detail below are summarized in Table ^. 
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Table 3: IIB Classification for QCD3 



To count the number of independent fluctuations for a field of given spin, we again imagine 
harmonic plane waves propagating in the AdS radial direction, r, with Euclidean time, 
X3. For example, the metric fluctuations in AdS^ 

G^u=g^u + h^y{x), (17) 
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in the fixed background g^y are taken to be of the form /i^jy(r, X3). There is no dependence 
on the spatial co-ordinates, xi = (xi,rE2), or the compactified "temperature" direction, r. 

3.1.1 Metric fluctuations 

A graviton has two polarization indices. If we were in flat space time, we could go to a 
gauge where these indices took values only among (xi, X2,t) and not from the set (r, 2:3). 
The polarization tensor should also be traceless. This leaves (3x4)/2 — 1 =5 independent 
components. In the AdS space time, we can count the number of graviton modes the same 
way, though the actual modes that we construct will have this form of polarization only 
at r ^ 00; for finite r, other components of the polarization will be constrained to acquire 
nonzero values ||lO[| . 

Therefore, a set of five independent polarization tensors can be characterized by the fol- 
lowing non- vanishing components at r —> 00. In the AdS^ black hole background, r is 
compact, so the rotations group is SO (2) in (xi, X2) and the five states are in 3 irreducible 
representations: A spin-2 doublet (helicities ±2), spanned by 

Gij : hi2 = h2i 7^ 0, and hu = -/i22 7^ , (18) 

a spin-1 doublet (helicities ±1), spanned by 

Gir ■ hrl = hir ^ 0, and hr2 = /i2,r 7^ (19) 

and a spin-0 state, 

Grr ■■ hrr = -2/lll = -2/l22 / . (20) 

These fluctuations are denoted by T3, V3, and 5*3 respectively in Table ^. 

3.1.2 Two-form fields 

Each 2-form field in AdS^ satisfies a field equation that includes a topological mass term. 
The 2-forms B^^ and C^y can be combined into one complex 2-form field B^y = B^^+iC^^. 
The field equation for B^y can be factorized into two first order equations, and each can 
be iterated leading to a second order equation of the form 

Max B^y + m\^sB^u = , 



13 



where Max is the Maxwell operator on 2-forms. For modes that have no dependence on 
the co-ordinates, one equation is massive, with m\^g = 16, and the other is massless. 
It can be shown that the massless equation has only pure gauge solutions; so they can be 
ignored. (See Appendix A for further details.) 

For the purpose of counting modes, polarizations for a massless 2-form in AdS^ can also be 
restricted to be transverse, with fields depending only on (x3,r). Therefore the polariza- 
tion tensor is an antisymmetric 2-tensor in the directions xi, X2,t, leading to 3 independent 
components. On the other hand, for a general massive 2-form field, longitudinal polariza- 
tions are allowed, so that the polarization is an antisymmetric tensor in the coordinates 
xi,X2,T,r, with 6 independent components. 

Nevertheless, the number of independent components for our massive 2-form B^i, is only 3 
(complex), as if we are dealing with a massless case. This is due to the fact that the second 
order equation above stems from a first order equation, e^u"^'' d^crB\p-\+ ^iB = 0, relating 
real and imaginary parts and leading to additional constraints. For instance, with fields 
depending on (x3,r), if we start with the polarization tensor having B12 = —B21 7^ , 
then the first order equation will constrain us to have specific values for Crr = —Crr 
and = —Crs 7^ 0, while allowing all other components of B and C to be zero. Thus we 
count as independent fields -B12, -Bit, -B2T, C12, Cit, C2t, i-e., there are three independent 
solutions for B^i, and three solutions for C^^. 

The non- vanishing polarizations of the B^u tensor can be grouped as: 

Bij-. Bi2 = -B2i / 0, (21) 

corresponding to spin-0 and 

Bir : Bir = —Brl 0, 

B2t = —Br2 / 0. 

corresponding to a spin-1 doublet. In Table ^, these fluctuations are denoted 
M3 respectively. They are degenerate with ones for the R-R 2-form C^^. 

3.1.3 Scalar fields 

In general, for each scalar field, there is a unique field equation with the plane wave 
dependence which we have been considering. There are three such scalar modes: The 
volume fluctuations in , (with m^^^ = 32), denoted by L3 and fluctuations for 



(22) 
by N3 and 
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the dilaton (j) and the axion C. However, as we have shown in an earUer paper |10|, the 
latter two spectra are degenerate with the 2"'""'" tensor fluctuations. Therefore, separate 
equations are not required 11, |l^] . 



3.2 Wave Equations and QCD^ Glueball Spectrum 



Metric fluctuations for QCD^ have been obtained previously by analyzing the linearized 
Einstein equations about the AdS^ x black hole background which leads to three 



independent equations, T3, V3 and 53 [10, 11, 12]. Fluctuations N^, M3 and L3 can be 
found similarly, leading to all together six independent equations for QCD-^. From the 
equation of motion, we determine the discrete spectrum with = im. (See Appendix A 
for details.) The full set of independent equations are: 
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1 



r)^3(r) 
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= 0, 
= 0, 



4^3 



1 



)M3(r) = 



(23) 



Each equation can be expressed in a variety of forms, depending on the choice of normal- 
ization. The following choices have been made so that each equation takes on a manifestly 
hermitian form: 

Tensor: 

hi, = qijr^T:,{r)e-^-\ (24) 
where i,j = 1,2, with qij an arbitrary constant traceless-symmetric 2x2 matrix. 



Vector: 



hir = qi^r^-W^{r)e-'^^^, 



(25) 



where qi is a constant 2-vector. 
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Equation: 


^3 


V3 


S3 


N3 


Ms 


L3 




2++/0 +/0++ 


t 1 1 

1++ 


^ 1 1 

0++ 


0+^/0 


1+ /I 


0++ 


n = 


11.588 


18.677 


5.457 


26.097 


44.272 


54.932 


n = 1 


34.527 


47.499 


30.442 


61.159 


84.318 


100.628 


n = 2 


68.975 


87.720 


65.123 


107.308 


135.932 


157.933 


n = 3 


114.910 


139.436 


111.141 


164.829 


199.073 


226.773 


n = 4 


172.331 


202.623 


168.601 


233.791 


273.717 


307.123 


n = 5 


241.237 


277.283 


237.528 


314.215 


359.861 


398.971 


n = 6 


321.627 


362.779 


317.931 


406.112 


457.507 


502.311 


n = 7 


413.501 


461.121 


409.815 


509.486 


566.631 


617.140 


n = 8 


516.860 


568.462 


513.180 


624.341 


687.279 


743.457 


n = 9 


631.703 


689.156 


628.028 


750.677 


819.329 


881.260 



Table 4: The mass spectrum, m^, for QCD3 Glueballs 



Scalar: This case has been treated carefuhy using several different gauge choices [10, 12]. 
Here we adopt the form suggested by Constable and Meyer |12| with 



hr 



(r^ — r ^)S'3(r)e 



(26) 



Two- forms and volume scalar: For B12, consider perturbations of the form 

Si2 = r2Ar3(r)e*'=3X3^ (27) 

with Bir = B2t = 0. Alternatively, we consider 

(28) 



3[r)e 

with B12 = 0, where qi is an arbitrary constant 2-vector. Lastly, for the volume perturba- 
tion, we consider 

h'^ = L3{r)e-^^\ 

To calculate the discrete spectrum for each of these equations one must again apply the 
correct boundary conditions at r = 1 and r = 00 as mentioned in the case of QCD4, earlier. 
(See Appendix B for a listing of boundary conditions.) The resultant spectrum is given 
in Table |. 



Similarly, we have also compared our numerical QCD3 masses to the WKB approxima- 
tions, 

~ ^2(^2 ^ ^ ^ ^ ^) ^ Q(l^ ^ (29) 

n 
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where = 167r(r(3/4)/r(l/4))^, with integer constants, 6, hsted in Table ^, determined 
analyticahy. Again, the constants 7 were fits to the numerical data. (See Appendix B). 



4 Parity and Charge Conjugation assignments 

Next we determine how the supergravity fields and therefore the glueballs couple to the 
boundary gauge theory. This allows us to unambiguously assign the correct parity and 
charge quantum numbers to the glueball states. For this purpose we consider the effective 
Born-Infeld action on the branes. 

4.1 QCD4 

The 4-d gauge theory for QCD4 is obtained by dimensional reduction from a 5-d gauge 
theory, which is the low energy dynamics of D4-branes in 10-dimensional Type IIA string 
theory. Although this 10-d theory may itself be regarded as a dimensional reduction of 
11-d M-theory for membranes, it is sufficient and more convenient to consider the 10-d 
theory itself. 

Since supergravity fields can be thought of as coupling constants for gauge theory op- 
erators, their quantum numbers can be assigned by the parity and charge conjugation 
invariance of the overall action, (supergravity field times composite operator). For sim- 
plicity let us consider the coupling of a supergravity field to just one D4-brane — this 
coupling is given by a Born-Infeld action plus a Wess-Zumino term, 

S = j d^xdet[Gf,^ + e-<^/2(^^^ ^ ^ J d*x{CiF A F + C3 A F + C5) , (30) 

where fj,,!^ = 1,2,3,4, r. Later we will argue that our quantum number assignment is 
correct also for the non-abelian case of N coincident D branes. In the 5-d field theory, 
we have the space-time world volume co-ordinate xi, X2, X3, X4, r with r compactified on 
S^. The Euclidean time coordinate we take to be X4. After dimensional reduction the 
physical fields will be characterized by their representation under the little group SO (3) 
of rotations on the spatial co-ordinates Xi, i = 1,2, 3, in the 4-d theory. 

For the 5-d gauge fields, we define parity by 

P : Ai[xi, X4^,t) > —Ai[—Xi,X4^,T), 



17 



P : A4{Xi,X4,T)^ A4{-Xi,X4,T), 

P : Ar{Xi,X4,T)^ Ar{-Xi,X4,T) , (31) 

for and T ^ T. For the Euclidean space, this is the only discrete 

symmetry. However after compactification to x S"*", we can define another parity (not 
related by a 5-d proper Lorentz transformation) by inverting the r co-ordinate on S"*". 
Thus we define a separate discrete r— parity transformation P,- : r — > — r, 



Pr 
Pr 
Pr 



Aj(Xj,X4,T) > Ai[xi,X4, t), 
Ai{xi,X4,T) A4{xi,X4,-T), 

Ar{Xi, X4, t) -AriXi, X4, -t) . (32) 



Charge conjugation for a non-abelian gluon field is 

C : \T,Al{x) ^ -kKAlix) (33) 

where are the Hermitian generators of the group. In terms of matrix fields {A = iTaA"), 
C : Aij,{x) — —A'j^{x). This leads to a subtlety. For example consider the transformation 
of a trilinear gauge invariant operators, 

C : Tr[F^-^^i^-^Fi^^i^^F^^^^] — rr[Py^3 1/3^^21^2 -^mi^^i] • i^^) 

The order of the fields is reversed. Hence the symmetric products, d^^'^Ff F^^l, have 
C = —1 and the antisymmetric products, /"^"^FfFl-^g , C = +1. Of course using a single 
branc, we can only find symmetric products. For reasons explained further in Sec. 5, we 
will only encounter symmetric traces over polynomials in F, designate by Sym Tr[F^,y • • •]. 
Even polynomials have C = +1 and odd polynomials C = — 1. 



4.1.1 Graviton couplings 

Expanding the Born-Infeld action, we can now read off the J^^(Pr) assignments: 
iFrom the coupling, G^^T'^'' - Gf,^Tr{F^xF^) + ■■■ , we obtain 

Gij ^ 2++ {Pr = +), Gir ^ 1-+ (Pr = +) G^r ^ 0++ {Pr = +) . (35) 

Under compactification of 11-d supergravity theory, G^,ii becomes the Ramond-Ramond 
1-form C^, which couples as ~ e^^'^'^^^C^ Sym T]:[Fy\F^j^W] , where W is an an even power 
of fields F. Consequently, the coupling, e^^^Ci Sym Th:[FrjF}^4W] + • • • leads to 

Ci ^ 1++ {Pr = -) . (36) 
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Similarly, e^^''CrTv{FijF4kW) + ■■■ gives 

a ^ 0-+ {Pr = +) , (37) 

and Gii^ii leads to the dilaton (p with coupling (pTiF"^ , 

4> ^ 0++ {Pr = +). (38) 

4.1.2 Two-form, three-form fields, and volume scalar 

Consider first the NS-NS 2- form field B^y . This field in the Type IIA theory arises from the 
3- form field of the 11-d supergravity theory when the components of the 3-form field are 
For U (1) gauge theory in leading order this field couples as B^^F^'^. More generally 
in the SU{N) gauge theory, we must have a multi-gluon coupling, Bj^n, SymTr[F^yW]. 
(Again W is an an even power of fields F and the trace is symmetrized.) 

To determine the parity, assume for the supergravity modes that we are in a gauge where 
the indices of the 2-form do not point along X4,xii,r. With i,j = 1,2,3, this leads to 
couphng Bij SymTr[F^^W] with 

Bij ^ 1 + ^ iPr = +), (39) 
and coupling Bi-r SymTr[F^'^W] with 

Bir ^ 1— {Pr = -) . (40) 

An analogous analysis can also be carried out for 3-form fields, Cijk and Cijr- The coupling, 
C7i23 SymTvlF^^W] leads to 

Cijk ^ 0+- {Pr = -) , (41) 
and coupling e^^^Cijr Symlh:[F^^W] leads to 

Cijr ^ 1— {Pr = +) . (42) 

Lastly, the volume scalar couples as /i" SymTrF'^ + • • • giving 

K ^ 0++ {Pr = +) . (43) 
The complete parity and charge conjugation assignments are given in Table 0. 
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4.2 QCDs 



The 3-d gauge theory is obtained by dimensional reduction from a 4-d gauge theory. To 
find the symmetries of the interactions, we consider the Born-Infeld action plus Wess- 
Zumino term, describing the coupling of a supergravity field to a single D3-brane, 

S = j (fx det[G^^ + e-^/2(^^^ ^ ^ j ^'^^(CoF A F + C2 A F + C^) , 

where = l,2,3,r. As in the case of QCD4, we will find the charge conjugation and 
parity assignments with the help of the symmetries of the 4-d gauge theory and then take 
the dimensional reduction to the 3-d theory after compactification of the coordinate r. 
The Euclidean time is taken to be X3 and the spatial co-ordinates, Xj, i = 1, 2. 

After dimensional reduction the physical fields will be characterized by their representation 
under the little group of the space co-ordinates of the 3-d theory: this is the group SO (2) 
rotations in the xi,X2 plane. However, unlike the case of QCD4, the usual spatial inversion, 
Xi — > —Xi with X3 — > X3 and r — > r, is a rotation R{9) in 5*0(2) with = vr; it therefore 
does not lead to a discrete symmetry. A discrete symmetry can be defined by xi — > xi and 
X2 —X2, as was pointed out in the lattice studies by Harte and Philipsen ||l^. However, 
the sole manifestation of this symmetry is the helicity "doublets", A = it J, for states with 
spin J > 0. We have already taken this degeneracy into account. 

On the other hand, r-parity remains a discrete symmetry of the action, as is the case for 
QCD4. In this paper, we shall define "parity" for QCD-^ as P = R{tt) x where Xi — > —Xi, 
X3 X3, and T —T. This of course is precisely the parity for the uncompactified 4-d 
theory with X3 treated as Euclidean time. 

4.2.1 Graviton, dilaton and axion states 

The graviton G^u couples as G^^T^'^ as in the case of QGD4. Because an even number of 
gluons occur in the field operators, the charge conjugation for all such states are C = +. 
For parity, we assume we are in a gauge where the indices of G^u do not point along X3, r. 
From the coupling, G^,y'Tt[F'^'^F^] + • • •, we get states 

GiJ ^ 2++ {Pr = +), G,r ^ 1++ {Pr = "), G^r ^ 0++ (P, = +) . (44) 

The dilaton couples as </>Tri^^, leading to 

^ 0++ (Pr = +), (45) 
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and for an axion coupling CoTr(Fi2-FT-3)) 

Co ^ 0-+ {Pr = -) . (46) 



4.2.2 Two-form fields and volume scalar 



Consider first the NS-NS 2- form field Bfj_i, . For parity, again assume that we are in a gauge 
where the indices of the 2-form do not point along X3,r. With i,j = 1,2, the coupling 
Bij SymTi[F'm] leads to 

Bij ^ 0+- {Pr = +) , (47) 

and Bir SymTr[F'''W] to 

Bir ^ 1+- {Pr = -). (48) 
Finally for the Ramond-Ramond 2-form Cfj_u, we have the coupling C12 SymTr[Fr3W], so 

C12 ^ 0^- {Pr = -) , (49) 

and e^^Cri SymTv[Fj-iW] giving 

C„ ^ 1~- {Pr = +)- (50) 

Finally, as in the case of QCD4, the volume scalar couples as /iJSJTrF^ + • • • so that 

^ 0++ {Pr = +) . (51) 

The complete parity and charge conjugation assignments are given in Table |3|. 



5 Discussion 



Lastly we turn to the question of how well the strong coupling limit for the Maldacena dual 
theory of QCD represents the infrared physics probed by the glueball spectra. Happily 
we now have a rather definitive lattice glueball spectrum by Morningstar and Pearson 
with which to make comparisons (See right side of Fig. ^ below). 
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5.1 Comparison with Lattice Glueball Spectrum 



Originally, claims were made about accurate comparisons to a few percent for isolated 
(scalar) mass ratios. As we pointed out in Ref. pO| ] for QCD^, the lowest mass scalar 
comes from the gravitational multiplet, not the dilaton. A similar spectrum is observed for 
QCD4. Consequently such accurate mass ratios were a misconception. This should not be 
regarded as a failure, since any reasonable expectation of a strong coupling approximation 
should not give quantitative results. On the other hand, there is a rather remarkable 
correspondence of the overall mass and spin structure between our strong coupling glueball 
spectrum and the lattice results at weak coupling for (see Fig. |2| below.) Apparently 

the spin structure of type IIA supergravity does resemble the low mass glueball spin 
splitting. The correspondence is sufficient to suggest that the Maldacena duality conjecture 
may well be correct and that further efforts to go beyond strong coupling are worthy of 
sustained eff'ort. 
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Figure 2: The AdS glueball spectrum for QCDj^ in strong coupling (left) compared with 
the lattice spectrum [^] for pure SU(3) QCD (right). The AdS cut-off scale is adjusted 
to set the lowest 2+"'" tensor state to the lattice results in units of the hadronic scale 
1/ro = 410 Mev. 
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Note that for each value of PC = (++,—+, H — , ), the lowest state is present in 

approximately the right mass range. In addition, the exact 2^+/0^+ degeneracy for 
AdS strong coupling corresponds to a relatively small splitting in the lattice calculations. 
Finally, there is a radial excitation of the pseudoscalar 0* that suggests that even this 
effect is approximated. This is intriguing because in the supergravity description the 
radial mode is a standing wave in the extra 5th dimension whereas in the lattice it is a 
conventional radial mode. Apparently scale changes in 4-d are being represented by the 
distance into the extra "warped" 5th axis. 



At higher masses the discrepancies increase. One reason is the obvious fact that on the 
supergravity side all orbital excitations of higher spin states are pushed to infinity in strong 
coupling by virtue of the divergent string tension. 

For example, the 3^^ state is a purely stringy effect outside of the classical limit of 
supergravity. 



Finally, we must emphasize that our comparison is premised on the neglect of many 
"spurious" states in the strong coupling limit that are in the wrong superselection sector 
to survive in the conjectured weak coupling limit of QCD. For example, all the Kaluza- 
Klein modes in the compact thermal manifold and the sphere (or S^) have masses 
at the cut-off scale. (The first mode on the thermal circle has a KK mass scale = 4 

in type IIB and m'j^j^ = 9 in type IIA in the units used in Table § and |2| respectively.) 
But these spurious KK modes all carry conserved U{1) or R charges that are absent in the 
target theory. We assume they will disappear in the continuum limit. A subtler situation 
occurs in the QCD4 example. Because normal modes in the extreme strong coupling 
limit do not distinguish between the compact 11th dimension and the spatial co-ordinates 
xi,X2,X3 on the brane, the spectrum actually has an exact S'0(4) symmetry. Thus there 
are additional states (listed in Table || but ignored in Fig. ^) exactly degenerate with 
the physically reasonable 2+^/0^^, ^, 1^ and 1 states. They are all odd under 
the discrete symmetry of reflecting the thermal circle (i.e. Pr = —) so they also lie 
in another superselection sector. A major challenge is to understand how this SO{4) 
symmetry is lifted and if the unwanted states remain at the cut-off in weak coupling. 
Physically the 11th axis is very different. The membranes of 11-d M-theory wraps this 
axis. Another possibility worth exploring is modifying the background metric with an 
orbifold that projects directly onto the even r— parity sector for QCD. 
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5.2 Constituent Gluon Picture 



The basic idea behind the AdS/CFT correspondence in the context of glueballs is similar 



to an observation made much earher by Fritzsch and Minkowski |22], by Bjorken |23] 
and by Jaffe, Johnson and Ryzak |2^. Namely that the low mass glueball spectrum can 
be qualitatively understood in terms of local gluon interpolating operators of minimal 
dimension. 

For example, Ref. |24] lists all gauge invariant operators for dimension d = 4,5 and 6. 
Eliminating operators that are zero by the classical equation of motion and states that 
decouple because of the conservation of the energy momentum tensor, the operators are in 
rough correspondence with all the low mass glueballs states, as computed in a constituent 
gluon [25 1 or bag model. Indeed more recently Kuti has pointed out that a more 



careful use of the spherical cavity approximation even gives a rather good quantitative 
match to the lowest 11 states in the lattice spectrum. 

Consequently it is interesting to compare this set of operators with the supergravity model. 
We list below all the operators for d < 6 except the operators with explicit derivatives 
(e.g. Tr[FDF]) and Tr[FDDF] ): 



Dimension 


State 


Operator 


Supergravity 


d = 4 


0++ 


Tr{FF) = E"" ■ E"" - B" ■ 




d = 4 


2++ 


Tij = Ef ■ Ef + Bf ■ Bf - trace 


Gij 


d = 4 


0-+ 


Tr{FF) = E"" ■B'' 


Cr 


d = 4 


0++ 


2 Too = E" ■ E" + B" ■ B" 


Gtt 


d = 4 


2-+ 


Ef ■ B^ + Bf ■ E^ - trace 


Absent 


d = 4 


2++ 


Ef ■ E^ - Bf • B^ - trace 


Absent 


d = 6 


(1,2,3)±- 


Tr{F^AFp^,Fxn}) ~ d^'^F^F^F'^ 


Bij , Cijr 


d = 6 


(0,1,2)±+ 


Tr{F^,[Fp,,Fx^]) ~ p^-F'^F^F^ 


Absent 



In this table we have used a Minkowski metric. The classification is parallel to our dis- 
cussion of couplings in the Born-Infeld action for D4 branes (see Sec. 4.1), except that all 
components and derivatives in the 5th, i.e., r, direction are zero and therefore all Pr = —1 
states are absent (See Table |l|). The column on the right lists the supergravity mode that 
couples (after dropping r components) to each operator. 

Several observations are in order. For the d = 4 operators, there is complete agreement 
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on the 3 lowest quantum numbers: J^^ = 0++,2++,0 ^. Also there is agreement on 
the absence of a low mass 1 state that Ref. [24| attributes to the conservation law that 
decouples the operator corresponding to the momentum tensor Tqi = E"" x B'^. In this 
context, it is worth commenting on the two independent sources for states — the 
condensate, Tr(FF), and the energy density. Too = ^(-E/" • E"" + B°' ■ B""). Naively one 
might conclude the second one coming from the conserved energy momentum tensor should 
be dropped, for the same reason we dropped the operator for momentum conservation, 
Toj. However in fact for QCD because of the conformal anomaly for the trace of the 
energy momentum tensor it is easy to show that the decoupling argument fails. A bag-like 
model circumvents the decoupling because the bag itself implies a scale breaking "vacuum" 
(empty bag) thus introducing an extra four- vector, the bag velocity u^. Our AdS black hole 
background, which is a key ingredient in our approach, also breaks conformal invariance. 
Consequently all agree that there is an extra low mass 0"''+ state in addition to the one 
which in our case is degenerate with the tensor 2++. The lattice data clearly favors this 
low mass 0^^ state, in agreement with our AdS spectrum. Finally, one rather low mass 
state, the 2 ^, is missing in the AdS spectrum. This state is clearly present in the lattice 
spectrum and is identified in the bag model. 



At d = 6, we have two states identified in the C = -1 symmetric trace, Tr{Ffj_y{Fpfj, Fxr^}): 
the l+~ state for operator d"-^'^Ba{Eb-Ec) and the 1 state for the operator d"-^'^Ea{Eh-Ec). 
These states are clearly related to the field content of a IIA supersymmetric multiplet. 
The higher spin representation are not present at strong coupling. Moreover, we have 
no states corresponding to the antisymmetric trace operators, Tr{F[F, F]). It appears 
that in the limit where we restrict to supergravity modes and ignore the massive stringy 
states, we will be able to obtain only the glueball state with the symmetric d"*'^ coupling 
between the group indices, and not the state with the antisymmetric /"^"^ coupling. If we 
consider the chiral primaries of the 4-d Yang-Mills theory (which was the theory on the 
boundary before we took the r direction to be compact), then we find that these had the 
form Sym Tt{X^X^ X'^) - i.e, we have a symmetric trace over the fields. Other operators 
that couple to the supergravity fields will be supersymmetry descendents of these chiral 
primaries, but the symmetry in the trace would be maintainecl|^. This fact may be related 
to the observation of Tseytlin |jl6| that generalizing the Born-Infeld action to a non-abelian 
case gives rise to symmetric trace operators in the field theory. 



One legitimate point of view is simply to suppose that all missing states must by definition 
be stringy effects that will be restored in weak coupling. However we prefer to look on 
this as a possible clue on constructing a better initial geometry for the supergravity/QCD 
duality proposals. 

'^We thank W. Taylor for a discussion on this point. 



25 



5.3 Strong coupling Expansion for Pomeron Intercept 



We shall end this discussion with a comment on the slope of the leading glueball trajectory 
as way to estimate the crossover value for the bare coupling, where continuum physics 
might begin to hold. The Pomeron is the leading Regge trajectory passing through the 
lightest glueball state with J^^ = 2^^. In a linear approximation, it can be parameterized 
by 

ap{t) = 2 + ap{t-ra\), (53) 
where we can use the strong coupling estimate for the lightest tensor massQ, 

1 



mT=^ [9.86 + 0(^)]/3-i . (54) 



Moreover if we make the standard assumption that the closed string tension is twice that 



between two static quark sources [19|, we also have a strong coupling expression for the 
Pomeron slope, 

' ' +0(-7^)]/?'. (55) 



Putting these together, we obtain a strong coupling expansion for the Pomeron intercept, 

47r 1 

ap(0):.2-0.66(-^)+0(^). (56) 



Turning this argument around, we can estimate a crossover value between the strong and 

weak coupling regimes by fixing ap(0) ~ 1.2 at its phenomenological value |17]. In fact this 

yields for QCD^ at = 3 a reasonable value for astrong = 5^/471 = 0.176 for the crossover. 

Much more experience with this new approach to strong coupling must be gained before 

such numerology can be taken seriously. However, similar crude argument have proven 

to be a useful guide in the crossover regime of lattice QCD. One might even follow the 

general strategy used in the lattice cut-off formulations. Postpone the difficult question 

of analytically solving the QCD string to find the true UV fixed point. Instead work at a 

fixed but physically reasonable cut-off scale (or bare coupling) to calculate the spectrum. 

If one is near enough to the fixed point, mass ratios should be reliable. After all, the real 

benefit of a weak/strong duality is to use each method in the domain where it provides the 

natural language. On the other hand, clearly from a fundamental point of view, finding 

analytical tools to understand the renormalized trajectory and prove asymptotic scaling 

within the context of the gauge invariant QCD string would also be a major achievement 

— an achievement that presumably would include a proof of confinement itself. 

^ We have adopted the normalization in the AdS-black hole metric to simplify the coefficients, e.g., for 
AdS^ , Qtt = — r~'*. This corresponds to fixing the "thermal-radius" R\ = 1/3 so that (3 = 2tvRi = 27r/3. 
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A Wave Equations 



In this appendix we outline the derivation of the wave equations that were used to find 
the energy levels in the supergravity theory. First we take the case of QCD^, for which 
we have the metric, 

ds^ = {r'-}^)dT^+r' ^ dx^^+{r^-\y^dr' + dnl, (A.l) 

^ 1=1,2,3 ^ 

where X3 is the Euclidean time direction. 

The simplest equation is the scalar wave equation for the dilaton and the axion. At the 
linear perturbation level, both satisfy 



,^,M = _L[^_^5M-^]^^ = 0. (A.2) 



-9 



We introduce a plane wave ansatz, 

^ = Ur)e"'^^^ , (A.3) 
with zero momentum and mass, m = ik^ providing the equation (8) for T3 in the text. 

Fluctuations in the volume of the sphere S^, which is a fiber at each point of space time, 
provides another scalar mode. This scalar has an AdS mass squared equal to 32. The field 
equation is 

^ ' '^<?^'^^/^],.-32</. = 0, (A.4) 



V-9 

L^{r)e^^^^^ , which reduces to the equation (8) for L3 in the text. 



Next let us address the case of the two- form fields, B^y and C^y, which at the linear level 
are conveniently combined into a single complex field, Buy = Bny + iCny. It was shown 



in [27| that this field satisfies the equation, 

(Max - k{k + 4))S^^ + 2ie^yP^''dpBxa = , (A.5) 
for A; = 0, 1, • • • harmonics on 5^. The Maxwell operator is defined by 

Max5^^ = H^,yx'^ . (A.6) 

in terms of the field strength, 

H^yx = d^Byx + dyBx^, + dxBpy . (A.7) 
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Since the second order differential operator factorizes into two first order operators, solu- 
tions fall into two classes, 

{2kI + i*D)Bl^J =0, {2{k + 4)I-i*D)Bj^J = 0, (A.8) 

where {*DB)fj,i, = Cui.f^'^dpBxa^ and I is the identity matrix. 

It is convenient to iterate these first order equations to get the second order equations, 

(Max - = o, (Max - {k + 4f)Bj^J = . (A.9) 

We are interested in fields with no dependence on the coordinates of the sphere, so we can 
take k = 0. It can be shown that the first class of solutions, Bj;^J , are pure gauge, so we 
are only interested in the second class, which has an effective mass squared of 16 for the 
field Bj^J. 

As explained in the text, one must of course check that solutions to the second order 

~ (2) 

equation for B^iJ , actually are valid solution to the original wave equation. This reduce 
the number of independent tensor fields, B^i,, from 6 to 3. For example with the ansatz, 

Bi2 = N^{r)r'^e'^^''^ (A. 10) 

the first order equations will determine and B^r once we have a solution of the second 
order equation for Bi2- This does not place any constraints on the solution for B12 itself. 
We have defined B\2 in terms of the normalized coefficient N^{r) to obtain a hermitian 
field equation for similar to our earlier scalar mode Ts. In a similar manner we can find 
the equation for Bit- We can solve the second order equation for this fiuctuation without 
constraint, and then the requirement arising from the associated first order equations 
determines corresponding values of ^23 and B2r- After adopting the ansatz indicated in 
the main text, we again obtain a hermitian equation for M^{r). 

The graviton perturbations arise from the Einstein action with a cosmological constant, 
expanded around the given background. We write G^i, = g^y + h^i,. the equation for the 
perturbation h^^u is 

- ^ V;A'^ - \hi,^,. + \hf,xy^ + + 4V = (A.ll) 

Near the boundary at r = 00, we can choose a gauge to make the perturbations transverse 
to r, 0:3 and traceless. It turns out that we can maintain this condition for all r for 
perturbations of the form hi2-, and for perturbations of the form h\r- In these cases the 
above equation for h^^ gives immediately the wave equations to be solved. But keeping in 
mind the decomposition in spin eigenstates in the X1—X2 plane, we also find that we have to 
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consider a spin-0 perturbation which at infinity has the form Ht-t, with hn = /i22 = — ^/irr- 
In this case we can choose a gauge to make /i3^ = 0, but for finite r we wih find in this 
gauge that hrr 7^ and also that the part transverse to r, is not traceless. Thus we 
have to keep hrT^hn = /122, /irr as independent coupled functions in the analysis. It was 



shown in [10| how these equations can be reduced to one effective equation which can then 



be solved in the same way as the equations for the other fields. A nicer choice of gauge 



was used in [12| which led to an equivalent but simpler equation. We will use the latter 
source for the equations, especially since the results there include all dimensions, and so 
can be used for the case of QCD4, as well. 

Now we turn to the case of AdS''' x S^, which is very similar. The metric is 

ds^ = {r^-\)dT^ + r^ ^ dx} + {r^--^r^dr^ + -dnl, (A.12) 

^ 1=1,2,3,4,11 

where we have 2:4 as the time direction and we note that the radius of the S'^ is half the 
curvature radius of the AdS space: this will affect the masses arising from the deformations 
of the sphere. 

There is a three-form field A^^x which behaves in a manner similar to the two-form 
field B^u discussed above. Its field equation can again be factorized into two first order 
equations, which we iterated to second order equations. One factor at A; = corresponds 
to pure gauge, while the other at A: = has the value m^^g = A{k + 3)^ = 36. For the 
scalar mode due to fluctuations of the volume of the S*^, we get a rn\^g = 4 x 18 = 72. 
For these considerations we arrive at the wave equations (22) given in the text. 

Finally for comparison with P. van Nieuwenhuizen Ref |^], we note that they have scaled 
the radius of the to unity, instead of |. Our choice was made to keep the radius of 
the AdS''^ equal to unity, to make the comparison between AdS^ and AdS''^ more natural. 
This change in metric scales the squared masses by a factor of 4 relative to Ref. |28]. 



B WKB and Variational Estimates 



First we change variable from r to a; = and express all twelve equations in the standard 
Sturm-Liouville form, 

^~'Ix^^^^^ ~^ w(a;)}</)„(x) = ml a{x)(j)n{x) , 
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where t(x), w{x), anda{x) are generalized "tension", "external force" , and "mass-density" 
respectively. In our case, we have for QCD^: Ti{x) = {x^ — x) and ai{x) = |; for QCD^: 
Ts{x) = {x^ — x) and (73 (x) = j. Force densities w{x) for all 12 cases are listed in Table ^. 
We shall use this as our starting point for carrying out variational and WKB analyses. 



B.l Variational Estimates for mh 



Solving for eigenstates, {(/>„(x)} and their corresponding eigenvalues, {m'^}, is equivalent 
to finding stationery points of the following functional, 

^,,,_ irdx[T{x)<f>'{xr+wixmxr] 

^f^] jrdxaix)Hx? ' ^""-'^ 

with 771^ — r[(^72]* particular, the square of the mass for the lowest state, ttiq — r[<^o]; 
is the absolute minimum of T[(f)]. 



To be properly defined as a Sturm-Liouville problem, it is necessary to impose boundary 
conditions: t(x)(/>(x)(/)'(x) — > 0, for x ^ 1 and for x ^ 00, in accord with the boundary 
conditions stated earlier for our numerical solutions. Explicit limiting behaviors for all 12 
cases are listed in Table ^. 



Given a trial wave function, Eq. B.l provides a variational upper bound for rriQ. As 

we have shown in Ref. |10|, accurate variational estimates for ground-state masses can 
be obtained with minimum efforts. Here, we shall only attempt to obtain simple esti- 
mates by choosing trial functions so that integrals in the equation above can be evaluated 
analytically. 



The simplest possible trial wave function for each case can be chosen as a product of t{x) 
and indicated in Table |5[ Indeed, our variational approach has served us well by 

providing a useful consistency check for our numerical efforts along the way. These are 
also summarized in Table |5|. 



B.2 WKB 



As explained in Ref. [p^ , we begin a WKB analysis by first bringing our differential equa- 
tions from the Sturm-Liouville form into a radial-Schroedinger form , 

(f 

( - + V{x; m^))^{x) = EiPix), (B.2) 
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Equation 


w{x) 


X — > 00 


X — > 1 


Trial Function 


Bound for m-g 


T4 





—3 
X ^ 


1 


—3 
X ^ 


24.0 


V4 


■gX [X — ij 




(X — i j ' 


/'^3 i\l/2^-9/2 
(X — i j ' X ' 


r .0 


c 


— lUoX [ox —I) 




i 


x~^ 


o.io 


AT 

iV4 


^(OX — l)X 


X ' 


1 
i 


^-9/2 
X ' 


00. U 


IVI4 


(ox — 2,)\x — ij 


^-9/2 
X ' 


yx — I) ' 


(X — i j ' X 


HQ Q 

yo.o 


T 

Li 


1 0^2 


™ — 6 
X 


i 


^ — 6 
X 


1 on 


T3 





X 


1 


X ^ 


12.0 


V3 


X-l(x2 - 


X-2 


(X- 1)1/2 


(x2 - 1)1/2^-3 


20.0 


S3 


-16x(3x2 - l)-2 


X~2 


1 


X~2 


5.55 


N3 


(3x2 _ -^^^-1 


X^2 


1 


X^2 


27.0 


Ms 


x(3x2 -2)(x2 - 




(X- 1)1/2 


(x2 - l)l/2x-4 


46.7 


L3 


8x 


x-4 


1 


—4 

X 


56.0 



Table 5: Variational Estimates for QCD^^ and QCD3 



by scaling ^p{x) = ^r(x)0(x). Eigenvalues nin are found by solving for zero-energy bound 
states from below, E ^ 0^ . The potential V{x; w?) is given by 

To obtain the desired WKB estimate, we simply need to evaluate in the large m? limit 



the following integral: (n + l/2)7r = J'^ dx\J —V{x;m?), i.e., we seek a WKB condition 
in the form, (n + ^)7r = sq^ti + si + O(^), where coefficients sq and si can be explicitly 
evaluated. As explained in Ref. [|lO[, the "effective potential", y(x,m2), is y(x;m2) = 
l/(x;m2) + l/(4(x- 1)2). 

To isolate the m2-dependence, let us write y(x,m2) = m2Vo(x) + Vi(x), where Vo(x) = 
-cT(x)/r(x), i.e., vj'^^(x) = -l/(4(x3 - 1)), and v!f'\x) = -l/(4x(x2 - 1)). For instance, 
one immediately finds that sq = dx^ —Vq{x). For QCD4 and QCD3, they are s'^^ = 
^r(i)/r(|), and s[,^^ = ^r(i)/r(f ). The remaining piece, Vi{x), is listed in Table |. 

To find coefficient si, we need to know the behavior of the ratio of Vo(x) to Vi{x) near 
X — > 1 and X 00. The dominant behavior of Vi{x), in these limits, can be characterized 
by two indices, 

T^i(x)~^^ + 0(^), for x^l, (B.4) 
Vi(x) ~ ^ + 0(^), for x^oo. (B.5) 
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Equation # 


Effective Potential Vi{x) 


61 


6r 


6 
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1,2 9 
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10 
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1 1 3 1 

iMJ /irv.^ f ^'2 T '\'A 
4i{X— L) 4iX X [X — i) 
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3 


2.01 
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1 1 3 
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1 
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4 
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1 1 3 1 16 
4{a;-l)^ ' Ax'^ x-^(x^-\y^ (x^-l)(3x^-l)^ 





2 


3 


-3.70 


N3 


1 1 3 1 1 3x-^-l 
4{a;-l)^ ' Ax^ x^{x^-\y^ ' x^(x^-l) 





4 


5 


4.71 


M3 


1 1 3 1 3a;+l 
4(x-l)^ ' Ax^ ' x^{x+\y^{x-\) 


1 


4 


6 
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L3 


1 1 3 1 1 8 
4(x-l)^ ' Ax^ x'-'ix^-iy-' ' x^-1 





6 


7 


9.50 



Table 6: WKB indices for QCD4 and QCD3 



One finds that si = -(f + 6r), and = {^)'^{n^ + Sn + 7} + 0(n"^) where 

5 = l + (5, + 5„ (B.6) 
as exhibited in Table ^. The values for 7 are found by a fit to the numerical values of m^. 
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